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^ ' Wc show that unless P=NP. there exists no polynomial time (or even pseudo-polynomial time) 

, algorithm that can decide whether a multivariate polynomial of degree four (or higher even 

Q>^ ' degree) is globally convex. This solves a problem that has been open since 1992 when N. Z. 

Shor asked for the complexity of deciding convexity for quartic polynomials. We also prove that 
1. s' ' deciding strict convexity, strong convexity, quasiconvexity, and pseudoconvexity of polynomials 

, of even degree four or higher is strongly NP-hard. By contrast, we show that quasiconvexity 

\my ' and pseudoconvexity of odd degree polynomials can be decided in polynomial time. 
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The role of convexity in modern day mathematical programming has proven to be remarkably 
fundamental, to the point that tractability of an optimization problem is nowadays assessed, more 
often than not, by whether or not the problem benefits from some sort of underlying convexity. In 
^3 '• the famous words of Rockafellar [39] : 

. "In fact the great watershed in optimization isn't between linearity and nonlinearity, 

I but convexity and nonconvexity." 

But how easy is it to distinguish between convexity and nonconvexity? Can we decide in an efficient 
manner if a given optimization problem is convex? 

A class of optimization problems that allow for a rigorous study of this question from a com- 
putational complexity viewpoint is the class of polynomial optimization problems. These are op- 
^ ' timization problems where the objective is given by a polynomial function and the feasible set is 

■ described by polynomial inequalities. Our research in this direction was motivated by a concrete 

question of N. Z. Shor that appeared as one of seven open problems in complexity theory for 
numerical optimization put together by Pardalos and Vavasis in 1992 [36]: 

"Given a degree-4 polynomial in n variables, what is the complexity of determining 
whether this polynomial describes a convex function?" 

As we will explain in more detail shortly, the reason why Shor's question is specifically about 
degree 4 polynomials is that deciding convexity of odd degree polynomials is trivial and deciding 
convexity of degree 2 (quadratic) polynomials can be reduced to the simple task of checking whether 
a constant matrix is positive semidefinite. So, the first interesting case really occurs for degree 4 
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(quartic) polynomials. Our main contribution in this paper (Theorem 12. II in Section [2. 3p is to show 
that deciding convexity of polynomials is strongly NP-hard already for polynomials of degree 4. 

The implication of NP-hardness of this problem is that unless P=NP, there exists no algorithm 
that can take as input the (rational) coefficients of a quartic polynomial, have running time bounded 
by a polynomial in the number of bits needed to represent the coefficients, and output correctly on 
every instance whether or not the polynomial is convex. Furthermore, the fact that our NP-hardness 
result is in the strong sense (as opposed to weakly NP-hard problems such as KNAPSACK) implies, 
roughly speaking, that the problem remains NP-hard even when the magnitude of the coefficients 
of the polynomial are restricted to be "small." For a strongly NP-hard problem, even a pseudo- 
polynomial time algorithm cannot exist unless P=NP. See [H] for precise definitions and more 
details. 

There are many areas of application where one would like to establish convexity of polynomials. 
Perhaps the simplest example is in global minimization of polynomials, where it could be very useful 
to decide first whether the polynomial to be optimized is convex. Once convexity is verified, then 
every local minimum is global and very basic techniques (e.g., gradient descent) can find a global 
minimum — a task that is in general NP-hard in the absence of convexity |37j . [33]. As another 
example, if we can certify that a homogeneous polynomial is convex, then we define a gauge (or 
Minkowski) norm based on its convex sublevel sets, which may be useful in many applications. In 
several other problems of practical relevance, we might not just be interested in checking whether a 
given polynomial is convex, but to parameterize a family of convex polynomials and perhaps search 
or optimize over them. For example we might be interested in approximating the convex envelope 
of a complicated nonconvex function with a convex polynomial, or in fitting a convex polynomial 
to a set of data points with minimum error [29]. Not surprisingly, if testing membership to the set 
of convex polynomials is hard, searching and optimizing over that set also turns out to be a hard 
problem. 

We also extend our hardness result to some variants of convexity, namely, the problems of 
deciding strict convexity, strong convexity, pseudoconvexity, and quasiconvexity of polynomials. 
Strict convexity is a property that is often useful to check because it guarantees uniqueness of the 
optimal solution in optimization problems. The notion of strong convexity is a common assumption 
in convergence analysis of many iterative Newton-type algorithms in optimization theory; see, 
e.g., [H Chaps. 9-11]. So, in order to ensure the theoretical convergence rates promised by many of 
these algorithms, one needs to first make sure that the objective function is strongly convex. The 
problem of checking quasiconvexity (convexity of sublevel sets) of polynomials also arises frequently 
in practice. For instance, if the feasible set of an optimization problem is defined by polynomial 
inequalities, by certifying quasiconvexity of the defining polynomials we can ensure that the feasible 
set is convex. In several statistics and clustering problems, we are interested in finding minimum 
volume convex sets that contain a set of data points in space. This problem can be tackled by 
searching over the set of quasiconvex polynomials [29]. In economics, quasiconcave functions are 
prevalent as desirable utility functions [27| , [4] . In control and systems theory, it is useful at times to 
search for quasiconvex Lyapunov functions whose convex sublevel sets contain relevant information 
about the trajectories of a dynamical system [10], [3]. Finally, the notion of pseudoconvexity is 
a natural generalization of convexity that inherits many of the attractive properties of convex 
functions. For example, every stationary point or every local minimum of a pseudoconvex function 
must be a global minimum. Because of these nice features, pseudoconvex programs have been 
studied extensively in nonlinear programming [30], |12] . 

As an outcome of close to a century of research in convex analysis, numerous necessary, 
sufficient, and exact conditions for convexity and all of its variants are available; see, e.g., [H 
Chap. 3], [32], [17], [l3], [27], [3l] and references therein for a by no means exhaustive list. Our 
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results suggest that none of the exact characterizations of these notions can be efficiently checked 
for polynomials. In fact, when turned upside down, many of these equivalent formulations reveal 
new NP-hard problems; see, e.g., Corollarv 12.61 and 12. 8i 

1.1 Related Literature 

There are several results in the literature on the complexity of various special cases of polynomial 
optimization problems. The interested reader can find many of these results in the edited volume 
of Pardalos [35] or in the survey papers of de Klerk |15j . and Blondel and Tsitsiklis [7j. A very 
general and fundamental concept in certifying feasibility of polynomial equations and inequalities 
is the Tarski-Seidenberg quantifier elimination theory [41], [40], from which it follows that all of the 
problems that we consider in this paper are algorithmically decidable. This means that there are 
algorithms that on all instances of our problems of interest halt in finite time and always output the 
correct yes-no answer. Unfortunately, algorithms based on quantifier elimination or similar decision 
algebra techniques have running times that are at least exponential in the number of variables |5] , 
and in practice can only solve problems with very few parameters. 

When we turn to the issue of polynomial time solvability, perhaps the most relevant result for 
our purposes is the NP-hardness of deciding nonnegativity of quartic polynomials and biquadratic 
forms (see Definition 12. 2p : the reduction that we give in this paper will in fact be from the latter 
problem. As we will see in Section 12. 3( it turns out that deciding convexity of quartic forms is 
equivalent to checking nonnegativity of a special class of biquadratic forms, which are themselves a 
special class of quartic forms. The NP-hardness of checking nonnegativity of quartic forms follows, 
e.g., as a direct consequence of NP-hardness of testing matrix copositivity, a result proven by Murty 
and Kabadi [33]. As for the hardness of checking nonnegativity of biquadratic forms, we know of 
two different proofs. The first one is due to Gurvits [21j . who proves that the entanglement problem 
in quantum mechanics (i.e., the problem of distinguishing separable quantum states from entangled 
ones) is NP-hard. A dual reformulation of this result shows directly that checking nonnegativity 
of biquadratic forms is NP-hard; see |16] . The second proof is due to Ling et al. [28], who use 
a theorem of Motzkin and Straus to give a very short and elegant reduction from the maximum 
independent set problem in graphs. 

The only work in the literature on the hardness of deciding polynomial convexity that we are 
aware of is the work of Guo on the complexity of deciding convexity of quartic polynomials over 
simplices |20j . Guo discusses some of the difficulties that arise from this problem, but he does not 
prove that deciding convexity of polynomials over simplices is NP-hard. Canny shows in [9] that 
the existential theory of the real numbers can be decided in PSPACE. From this, it follows that 
testing several properties of polynomials, including nonnegativity and convexity, can be done in 
polynomial space. In [34], Nie proves that the related notion of matrix convexity is NP-hard for 
polynomial matrices whose entries are quadratic forms. 

On the algorithmic side, several techniques have been proposed both for testing convexity of sets 
and convexity of functions. Rademacher and Vempala present and analyze randomized algorithms 
for testing the relaxed notion of approximate convexity [38] . In [26] , Lasserre proposes a semidefinite 
programming hierarchy for testing convexity of basic closed semialgebraic sets; a problem that we 
also prove to be NP-hard (see Corollarv 12. 8p . As for testing convexity of functions, an approach 
that some convex optimization parsers (e.g., CVX [19]) take is to start with some ground set of 
convex functions and then check whether the desired function can be obtained by applying a set 
of convexity preserving operations to the functions in the ground set |14] . [8l p. 79]. Techniques 
of this type that are based on the calculus of convex functions are successful for a large range 
of applications. However, when applied to general polynomial functions, they can only detect a 
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subclass of convex polynomials. 

Related to convexity of polynomials, a concept that has attracted recent attention is the alge- 
braic notion of sos- convexity (see Definition I2.4p [22], [21], [25], [3], [29], [10]. This is a powerful 
sufficient condition for convexity that relies on an appropriately defined sum of squares decomposi- 
tion of the Hessian matrix, and can be efficiently checked by solving a single semidefinite program. 
However, in [T], (2], Ahmadi and Parrilo gave an explicit counterexample to show that not every 
convex polynomial is sos-convex. The NP-hardness result in this work certainly justifies the exis- 
tence of such a counterexample and more generally suggests that any polynomial time algorithm 
attempted for checking polynomial convexity is doomed to fail on some hard instances. 

1.2 Contributions and organization of the paper 

The main contribution of this paper is to establish the computational complexity of deciding con- 
vexity, strict convexity, strong convexity, pseudoconvexity, and quasiconvexity of polynomials for 
any given degree. (See Table 1 in Section [5] for a quick summary.) The results are mainly divided 
in three sections, with Section [2] covering convexity. Section [3] covering strict and strong convexity, 
and Section U] covering quasiconvexity and pseudoconvexity. These three sections follow a similar 
pattern and are each divided into three parts: first, the definitions and basics, second, the degrees 
for which the questions can be answered in polynomial time, and third, the degrees for which the 
questions are NP-hard. 

Our main reduction, which establishes NP-hardness of checking convexity of quartic forms, is 
given in Section 12.31 This hardness result is extended to strict and strong convexity in Section 13. 3^ 
and to quasiconvexity and pseudoconvexity in Section 14.31 By contrast, we show in Section 14.21 
that quasiconvexity and pseudoconvexity of odd degree polynomials can be decided in polynomial 
time. Finally, a summary of our results and some concluding remarks are presented in Section [5j 

2 Complexity of deciding convexity 



2.1 Definitions and basics 

A (multivariate) polynomial p{x) in variables x := 
a finite linear combination of monomials: 



(xi, . . . , Xn)"^ is a function from M" to M that is 



Pj^) — ^ '] C-gX — 'y ^ Cai,...,a„X]^ ^ • • • Xjj" , (1) 
a ai,...,an 

where the sum is over n-tuples of nonnegative integers Oj. An algorithm for testing some property 
of polynomials will have as its input an ordered list of the coefficients Cq. Since our complexity 
results are based on models of digital computation, where the input must be represented by a finite 
number of bits, the coefficients Cq, for us will always be rational numbers, which upon clearing 
the denominators can be taken to be integers. So, for the remainder of the paper, even when not 
explicitly stated, we will always have Cq € Z. 

The degree of a monomial is equal to ai -|- • • • + a„. The degree of a polynomial p{x) is 
defined to be the highest degree of its component monomials. A simple counting argument shows 
that a polynomial of degree d in n variables has ("'^'^) coefficients. A homogeneous polynomial (or 
a form) is a polynomial where all the monomials have the same degree. A form p{x) of degree d is 
a homogeneous function of degree d (since it satisfies p{\x) = A'^p(x)), and has i^'^'jf^) coefficients. 

A polynomial p{x) is said to be nonnegative or positive semidefinite (psd) if p{x) > for all 
X £ R". Clearly, a necessary condition for a polynomial to be psd is for its degree to be even. 
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We say that p{x) is a sum of squares (sos), if there exist polynomials qi{x), . . . , qm{x) such that 
p{x) = ^™iQ'f(a^)- Every sos polynomial is obviously psd. A polynomial matrix P{x) is a matrix 
with polynomial entries. We say that a polynomial matrix P[x) is PSD (denoted P{x) ^ 0) if it is 
positive semidefinite in the matrix sense for every value of the indeterminates x. (Note the upper 
case convention for matrices.) It is easy to see that P{x) is PSD if and only if the scalar polynomial 
P{x)y in variables (x; y) is psd. 

We recall that a polynomial p{x) is convex if and only if its Hessian matrix, which will be 
generally denoted by H{x), is PSD. 

2.2 Degrees that are easy 

The question of deciding convexity is trivial for odd degree polynomials. Indeed, it is easy to check 
that linear polynomials {d = 1) are always convex and that polynomials of odd degree d > 3 can 
never be convex. The case of quadratic polynomials {d = 2) is also straightforward. A quadratic 
polynomial p{x) = ^x'^Qx + q^x + c is convex if and only if the constant matrix Q is positive 
semidefinite. This can be decided in polynomial time for example by performing Gaussian pivot 
steps along the main diagonal of Q [33] or by computing the characteristic polynomial of Q exactly 
and then checking that the signs of its coefficients alternate |23| p. 403]. 

Unfortunately, the results that come next suggest that the case of quadratic polynomials is 
essentially the only nontrivial case where convexity can be efficiently decided. 

2.3 Degrees that are hard 

The main hardness result of the paper is the following theorem. 

Theorem 2.1. Deciding convexity of degree four polynomials is strongly NP-hard. This is true 
even when the polynomials are restricted to be homogeneous. 

We will give a reduction from the problem of deciding nonnegativity of biquadratic forms. We 
start by recalling some basic facts about biquadratic forms and sketching the idea of the proof. 

Definition 2.2. A biquadratic form b{x;y) is a form in the variables x = (xi,...,x„)^ and 
y = (yi, . . . , ym)^ that can be written as 



Note that for fixed x, b{x; y) becomes a quadratic form in y, and for fixed y, it becomes a 
quadratic form in x. Every biquadratic form is a quartic form, but the converse is of course not 
true. It follows from a result of Ling et al. [28] that deciding nonnegativity of biquadratic forms is 
strongly NP-hard even when n = m. In the sequel, we will always have n = m. 

It is not difficult to see that any biquadratic form b{x; y) can be written in the form 



(or of course as x B{y)x) for some symmetric polynomial matrix A{x) whose entries are quadratic 
forms. Therefore, it is strongly NP-hard to decide whether a symmetric polynomial matrix with 
quadratic form entries is PSD. One might hope that this would lead to a quick proof of NP-hardness 
of testing convexity of quartic forms, because the Hessian of a quartic form is exactly a symmetric 
polynomial matrix with quadratic form entries. However, the major problem that stands in the 




(2) 



i<j, k<l 



(3) 
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way is that not every polynomial matrix is a valid Hessian. Indeed, if any of the partial derivatives 
between the entries of A{x) do not commute (e.g., if ^"^g^^^^ 7^ ^^dxi^^ ) ' then A{x) cannot be the 
matrix of second derivatives of some polynomial. This is because all mixed third partial derivatives 
of polynomials must commute. 

Our task is therefore to prove that even with these additional constraints on the entries of A{x), 
the problem of deciding positive semidefiniteness of such matrices remains NP-hard. We will show 
that any given symmetric n x n matrix A{x), whose entries are quadratic forms, can be embedded 
in a 2n X 2n polynomial matrix H(x, y), again with quadratic form entries, so that H{x, y) is a valid 
Hessian and A{x) is PSD if and only if H(x, y) is. In fact, we will directly construct the polynomial 
f{x,y) whose Hessian is the matrix H(x,y). This is done in the next theorem, which establishes 
the correctness of our main reduction. Once this theorem is proven, the proof of Theorem 12.11 will 
become immediate. 

Theorem 2.3. Given a biquadratic form b{x;y), define the the n x n polynomial matrix C{x,y) 
by setting 

and let 7 be the largest coefficient, in absolute value, of any monomial present in some entry of the 
matrix C{x,y). Let f be the form given by 

2 n n 

f{x,y):=h{x-y) + '^(j2< + T.yt+ E + E vht)- (5) 

i=l i=l i,j=l,...,n i,j=l,...,n 

i<j i<j 

Then, b{x; y) is psd if and only if f{x, y) is convex. 

Proof. Before we prove the claim, let us make a few observations and try to shed light on the 
intuition behind this construction. We will use H{x, y) to denote the Hessian of /. This is a 
2n X 2n polynomial matrix whose entries are quadratic forms. The polynomial / is convex if and 
only if H{x, y)z is psd. For bookkeeping purposes, let us split the variables 
where Zx and Zy each belong to M". It will also be helpful to give a name to the second group of 
terms in the definition of f{x, y) in ([5]). So, let 

2 n. n 

9i-,y)--=^{E-' + Eyt+ E E y^yfj- (6) 

i=l 2=1 i,j=l,...,n i,j=l,...,n 

i<j i<j 

We denote the Hessian matrices of b{x, y) and g{x, y) with Hb{x, y) and Hg{x, y) respectively. Thus, 
H{x,y) = Hi,{x,y) + Hg{x,y). Let us first focus on the structure of Hi,{x,y). Observe that if we 
define 

oyidyj 

^y'^A{x)y = b{x;y). (7) 



then A{x) depends only on x, and 
Similarly, if we let 
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then B{y) depends only on y, and 



1 



B{y)x = b{x; y). 



(8) 



From Eq. ([7]), we have that b{x;y) is psd if and only if A{x) is PSD; from Eq. ([8]), we see that 
b{x;y) is psd if and only if B{y) is PSD. 
Putting the blocks together, we have 



Hbix,y) 



B{y) C{x,y) 
C^{x,y) A{x) 



(9) 



The matrix C{x,y) is not in general symmetric. The entries of C{x,y) consist of square-free 
monomials that are each a multiple of Xiyj for some i, j, with 1 < i, j < n; (see ([2]) and ([1])). 
The Hessian Hg(x, y) of the polynomial g{x, y) in ([6]) is given by 



Hg{x,y) 



2 

n 7 



Hl\x) 






Hf{y) 



(10) 



where 



Hl\x) 



12xf + 2 ^ 

i=l 
4xiX2 



4X1X2 



i=l,...,n 



12x2 + 2 ^ 



4X2X„ 



1=1, ...,n 



4:Xn-lXn 12x2+2 ^ X^ 



i=l,...,n 



(11) 



and 



12y2 + 2 



42/12/2 



i^i 



4yiy2 



4yi2/n 



122/2 + 2 y2 



4yi2/n 

42/22/^ 



i=l,...,n 



42/„_i2/n 122/2+2 5^ yf 



i=l,...- 



(12) 



Note that all diagonal elements of Hg^{x) and Hg'^{y) contain the square of every variable xi, . . . , x„ 
and 2/1, • • • , 2/n respectively. 

We fist give an intuitive summary of the rest of the proof. If b{x;y) is not psd, then B{y) and 
A{x) are not PSD and hence Hh{x,y) is not PSD. Moreover, adding Hg{x,y) to Hh{x,y) cannot 
help make H{x,y) PSD because the dependence of the diagonal blocks of Hb{x,y) and Hg{x,y) 
on X and y runs backwards. On the other hand, if b{x;y) is psd, then Hi,{x,y) will have PSD 
diagonal blocks. In principle, Hi)(x,y) might still not be PSD because of the off-diagonal block 
C{x,y). However, the squares in the diagonal elements of Hg{x,y) will be shown to dominate the 
monomials of C{x,y) and make H(x,y) PSD. 
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Let us now prove the theorem formally. One direction is easy: if b{x; y) is not psd, then /(x, y) 
is not convex. Indeed, if there exist x and y in such that y) < 0, then 



y'^A{x)y = 2b{x;y)<0. 



Zx=0,x=x,y=0,Zy=y 

For the converse, suppose that b{x; y) is psd; we will prove that H{x, y)z is psd and hence 



f{x,y) is convex. We have 



z'H{x,y)z = z'yA{x)zy + ziBiy)z, + 2ziCix,y)zy + + Hf (y) Zy . (13) 

Because ZyA{x)zy and z'^B{y)zx are psd by assumption (see d?]) and ([H])), it suffices to show that 
z'^H{x,y)z — ZyA{x)zy — z'^B{y)zx is psd. In fact, we will show that z'^ H{x,y)z — ZyA{x)zy — 
z'^B{y)zx is a sum of squares. 

After some regrouping of terms we can write 



^ 7 Tt722/ 



z^H{x,y)z - ZyA{x)zy - z^B{y)zx =pi{x,y,z) +P2{x,Zx) +p3{y,Zy), 



where 



pi{x,y,z) = 2zJC(x,y)2y + n%(^ 



^x,i 



P2{x,Zx)= r?jz. 



T 



and 



i=l 

bx\ 

2xiX2 
2x\Xyi 

22/1 yn 



j=i 



2xiX2 

5Xn 



2Xn^—\Xj; 



22/1^2 



5yi 



2X2Xn 



2yiyn 
2y2yn 



(14) 
(15) 

(16) 



2yn-iyn 5yl 



'-V 



(17) 



We show that (jl4p is sos by showing that pi, p2) and p3 are each individually sos. To see that 
P2 is SOS, simply note that we can rewrite it as 



P2{x,Zx) = n^7 



. k=l k=l 

The argument for p3 is of course identical. To show that pi is sos, we argue as follows. If we 
multiply out the first term 2ZxC{x,y)zy, we obtain a polynomial with monomials of the form 

±'^/3i,j,k,iZx,kXiyjZy^i, (18) 

where < /3ij^k,i < 7) by the definition of 7. Since 

± 2/3ij^k,iZx,kXiyjZy^i + fiij^k.izl^k^'i + f3ij,k,iy'jZy i = /3ij^k,i{^x,kXi ± yjZy^if, (19) 

by pairing up the terms of 2z'^C{x,y)zy with fractions of the squared terms z^ ^x? and Zy^yj, we 
get a sum of squares. Observe that there are more than enough squares for each monomial of 
2z'^C{x,y)zy because each such monomial ±2j3i^j^k,iZx,kXiyjZy^i occurs at most once, so that each 
of the terms z"^ ^x? and z'^ ly^ will be needed at most times, each time with a coefficient of at 
most 7. Therefore, pi is sos, and this completes the proof. □ 
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We can now complete the proof of strong NP-hardness of deciding convexity of quartic forms. 

Proof of Theorem \2.1\ As we remarked earlier, deciding nonnegativity of biquadratic forms is 
known to be strongly NP-hard [28]. Given such a biquadratic form b{x;y), we can construct 
the polynomial f{x,y) as in Note that f{x,y) has degree four and is homogeneous. More- 
over, the reduction from b{x;y) to f{x,y) runs in polynomial time as we are only adding to b{x;y) 
2n + 2(2) new monomials with coefficient and the size of 7 is by definition only polynomially 
larger than the size of any coefficient of b{x; y). Since by Theorem 12.31 convexity of /(x, y) is equiv- 
alent to nonnegativity of b{x;y), we conclude that deciding convexity of quartic forms is strongly 



An algebraic version of the reduction. Before we proceed further with our results, we make 
a slight detour and present an algebraic analogue of this reduction, which relates sum of squares 
biquadratic forms to sos-convex polynomials. Both of these concepts are well-studied in the liter- 
ature, in particular in regards to their connection to semidefinite programming; see, e.g., [28], [1], 
and references therein. 

Definition 2.4. A polynomial p{x), with its Hessian denoted by H{x), is sos-convex if the poly- 
nomial y'^H{x)y is a sum of squares in variables (x;y)^ 

Theorem 2.5. Given a biquadratic form b{x;y), let f{x,y) be the quartic form defined as in 
Then b{x;y) is a sum of squares if and only if f{x,y) is sos-convex. 

Proof. The proof is very similar to the proof of Theorem 12.31 and is left to the reader. □ 

Perhaps of independent interest, Theorems I2.3l and l2.5l imply that our reduction gives an explicit 
way of constructing convex but not sos-convex quartic forms (see [1]), starting from any example 
of a psd but not sos biquadratic form (see [TT]). 

Some NP-hardness results, obtained as corollaries. NP-hardness of checking convexity of 
quartic forms directly establishes NP-hardnes^ of several problems of interest. Here, we mention 
a few examples. 

Corollary 2.6. It is NP-hard to decide nonnegativity of a homogeneous polynomial q of degree 
four, of the form 



for some homogeneous quartic polynomial p. 

Proof. Nonnegativity of q is equivalent to convexity of p, and the result follows directly from 



NP-hard. 



□ 




Theorem 12. 1[ 



□ 



Definition 2.7. A set S C is basic closed semialgebraic if it can be written as 



S = {x £ M"| fi{x) > 0, i 



= l,...,m} 



(20) 



for some positive integer m and some polynomials fi{x). 



^See [3] for three other equivalent definitions of sos-convexity. 

■^All of our NP-hardness results in this paper are in the strong sense. For the sake of brevity, from now on we refer 
to strongly NP-hard problems simply as NP-hard problems. 
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Corollary 2.8. Given a basic closed semialgebraic set S as in i20\) . where at least one of the 
defining polynomials fi{x) has degree four, it is NP-hard to decide whether S is a convex set. 

Proof. Given a quartic polynomial p{x), consider the basic closed semialgebraic set 

£p = {{x,t) G]R"+i| t-p{x) >0}, 

describing the epigraph of p{x). Since p{x) is convex if and only if its epigraph is a convex set, the 
result follows H □ 

Convexity of polynomials of even degree larger than four. We end this section by extend- 
ing our hardness result to polynomials of higher degree. 

Corollary 2.9. It is NP-hard to check convexity of polynomials of any fixed even degree d > 4. 

Proof. We have already established the result for polynomials of degree four. Given such a degree 
four polynomial p{x) := p{xi, . . . , x„) and an even degree d > 6, consider the polynomial 

q{x,Xn+i) =p{x) 

in n + 1 variables. It is clear (e.g., from the block diagonal structure of the Hessian of q) that p{x) 
is convex if and only if q{x) is convex. The result follows. □ 

3 Complexity of deciding strict convexity and strong convexity 

3.1 Definitions and basics 

Definition 3.1. A function f : M" —?■ M is strictly convex if for all x ^ y and all A € (0, 1), we 
have 

f{Xx + (1 - A)y) < Xf{x) + (1 - X)f{y). (21) 

Definition 3.2. A twice differentiahle function f : R" — > M is strongly convex if its Hessian H{x) 
satisfies 

H{x) >z ml, (22) 

for a scalar m > and for all x. 

We have the standard implications 

strong convexity =^ strict convexity =^ convexity, (23) 
but none of the converse implications is true. 

3.2 Degrees that are easy 

From the implications in ()23p and our previous discussion, it is clear that odd degree polynomials 
can never be strictly convex or strongly convex. We cover the case of quadratic polynomials in the 
following straightforward proposition. 

Proposition 3.3. For a quadratic polynomial p{x) = ^x'^Qx + q'^x + c, the notions of strict 
convexity and strong convexity are equivalent, and can be decided in polynomial time. 

^Another proof of this corollary is given by the NP-hardness of checking convexity of sublevel sets of quartic 
polynomials (Theorem 14. 101 in Section [4. 3p . 
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Proof. Strong convexity always implies strict convexity. For the reverse direction, assume that p{x) 
is not strongly convex. In view of (I22p . this means that the matrix Q is not positive definite. If Q 
has a negative eigenvalue, p{x) is not convex, let alone strictly convex. If Q has a zero eigenvalue, 
let a; 7^ be the corresponding eigenvector. Then p{x) restricted to the line from the origin to x is 
linear and hence not strictly convex. 

To see that these properties can be checked in polynomial time, note that p{x) is strongly 
convex if and only if the symmetric matrix Q is positive definite. By Sylvester's criterion, positive 
definiteness of an n x n symmetric matrix is equivalent to positivity of its n leading principal minors, 
each of which can be computed in polynomial time. □ 

3.3 Degrees that are hard 

With little effort, we can extend our NP-hardness result in the previous section to address strict 
convexity and strong convexity. 

Proposition 3.4. It is NP-hard to decide strong convexity of polynomials of any fixed even degree 
d>A. 

Proof. We give a reduction from the problem of deciding convexity of quartic forms. Given a 
homogenous quartic polynomial p{x) := p{xi, . . . ,Xn) and an even degree d > 4, consider the 
polynomial 

q{x, Xn+i) := p{x) + x^,+i + ^{xj + • • • + 4 + xl_^^) (24) 

in n + 1 variables. We claim that p is convex if and only if q is strongly convex. Indeed, if p{x) is 
convex, then so is p{x) + x^_^^. Therefore, the Hessian of p{x) + x'^_^_i is PSD. On the other hand, 
the Hessian of the term ^(xf + • • • + + x^^^) is the identity matrix. So, the minimum eigenvalue 
of the Hessian of q{x,Xn+i) is positive and bounded below by one. Hence, q is strongly convex. 

Now suppose that p{x) is not convex. Let us denote the Hessians of p and q respectively by Hp 
and Hq. If p is not convex, then there exists a point x S M"" such that 

>^mm{Hp{x)) < 0, 

where Amin here denotes the minimum eigenvalue. Because p{x) is homogenous of degree four, we 
have 

Amm(-ffp(cx)) = C^\min{Hp{x)), 

for any scalar c S M. Pick c large enough such that X^i^{Hp{cx)) < 1. Then it is easy to see that 
Hq{cx, 0) has a negative eigenvalue and hence q is not convex, let alone strongly convex. □ 

Remark 3.1. It is worth noting that homogeneous polynomials of degree d > 2 can never be strongly 
convex (because their Hessians vanish at the origin) . Not surprisingly, the polynomial q in the proof 
of Proposition 13.41 is not homogeneous. 

Proposition 3.5. It is NP-hard to decide strict convexity of polynomials of any fixed even degree 
d>A. 

Proof. The proof is almost identical to the proof of Proposition 13.41 Let q be defined as in ()24p . 
If p is convex, then we established that q is strongly convex and hence also strictly convex. If p is 
not convex, we showed that q is not convex and hence also not strictly convex. □ 
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4 Complexity of deciding quasiconvexity and pseudoconvexity 



4.1 Definitions and basics 

Definition 4.1. A function f : M" — t- M is quasiconvex if its suhlevel sets 

S{a) := {x G M" I f{x) < a}, (25) 

for all a € M, are convex. 

Definition 4.2. A differentiable function f : M" — )> R is pseudoconvex if the implication 

Vf{xf{y-x)>0^f{y)>f{x) (26) 

holds for all x and y in M". 

The following implications are well-known (see e.g. [6l p. 143]): 

convexity =^ pseudoconvexity quasiconvexity, (27) 
but the converse of neither implication is true in general. 

4.2 Degrees that are easy 

As we remarked earlier, linear polynomials are always convex and hence also pseudoconvex and 
quasiconvex. Unlike convexity, however, it is possible for polynomials of odd degree d > 3 to be 
pseudoconvex or quasiconvex. We will show in this section that somewhat surprisingly, quasicon- 
vexity and pseudoconvexity of polynomials of any fixed odd degree can be decided in polynomial 
time. Before we present these results, we will cover the easy case of quadratic polynomials. 

Proposition 4.3. For a quadratic polynomial p{x) = ^x'^Qx + q^x + c, the notions of convexity, 
pseudoconvexity, and quasiconvexity are equivalent, and can be decided in polynomial time. 

Proof. We argue that the quadratic polynomial p{x) is convex if and only if it is quasiconvex. 
Indeed, if p{x) is not convex, then Q has a negative eigenvalue; letting x be a corresponding 
eigenvector, we have that p{tx) is a quadratic polynomial in t, with negative leading coefficient, so 
p{tx) is not quasiconvex, as a function of t. This, however, implies that p{x) is not quasiconvex. 

We have already argued in Section 12.21 that convexity of quadratic polynomials can be decided 
in polynomial time. □ 

4.2.1 Quasiconvexity of polynomials of odd degree 

In this subsection, we provide a polynomial time algorithm for checking whether an odd-degree 
polynomial is quasiconvex. Towards this goal, we will first show that quasiconvex polynomials of 
odd degree have a very particular structure (Proposition 14. 6p . 

Our first lemma concerns quasiconvex polynomials of odd degree in one variable. A version of 
this lemma is provided in [HI p. 99], without proof. We provide here a proof, for completeness. 

Lemma 4.4. Suppose that p{t) is a quasiconvex univariate polynomial of odd degree. Then, p{t) 
is monotonic. 
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Proof. Suppose that monotonicity does not hold for the points a < b < c. This can happen in two 
ways: 

p{a) < p{h) and p{c)<p{b), 

or 

p{a) > p{b) and p{c) > p{b). 

In the first case, pick a such that p{b) > a > max(p(a),p(c)). The sublevel set {x \ p{x) < a} 
includes a and c, but not b, and is therefore non-convex, which contradicts the quasiconvexity of p. 

In the second case, pick a such that min(p(a),p(c)) > a > p{b), and consider the sublevel set 
S{a) = {x I p{x) < a}. If the leading coefficient of p is positive, then S{a) contains contains b 
and some point to the left of a (since limt_^_ooP(^) = — oo) but not a, which is a contradiction. 
Similarly, if the leading coefficient of p is negative, then S{a) contains b and some point to the 
right of c but not c, which is a contradiction. □ 

Next, we use the preceding lemma to characterize the complements of sublevel sets of quasicon- 
vex polynomials of odd degree. 

Lemma 4.5. Suppose that p{x) is a quasiconvex polynomial of odd degree d. Then the set {x \ 
p{x) > a} is convex. 

Proof. Suppose not. In that case, there exist x,y,z such that z is on the line segment connecting 
X and y, and such that p{x),p{y) > a but p{z) < a. Consider the polynomial 

q{t) =p{x + t{y - x)). 

This is, of course, a quasiconvex polynomial with q{Q) = p{x), q{l) = p{y), and q{t') = p{z), for 
some t' G (0, 1). If q{t) has degree d, then, by Lemma [4.41 it must be monotonic, which immediately 
provides a contradiction. 

Suppose now that q{t) has degree less than d. Let us attempt to perturb x to x + x' , and y to 
y + y' , so that the new polynomial 

q{t) = p (x + x' + t{y + y' — X — x')) 

has the following two properties: (i) q{t) is a polynomial of degree d, and (ii) q{0) > q{t'), q{l) > 
q{t'). If such perturbation vectors x',y' can be found, then we obtain a contradiction as in the 
previous paragraph. 

To satisfy condition (ii), it suffices (by continuity) to take x',y' with small enough. 

Thus, we only need to argue that we can find arbitrarily small x',y' that satisfy condition (i). 
Observe that the coefficient of in the polynomial q(t) is a nonzero polynomial in x + x',y-\- y'; let 
us denote that coefficient as r(x + x',y + y'). Since r is a nonzero polynomial, it cannot vanish at 
all points of any given ball. Therefore, even when considering a small ball around (x, y) (to satisfy 
condition (ii)), we can find (x + x', y + y') in that ball, with r(x + x', y + y') ^ 0, thus establishing 
that the degree of q is indeed d. This completes the proof. 

□ 

We now proceed to a characterization of quasiconvex polynomials of odd degree. 

Proposition 4.6. Let p[x) be a polynomial of odd degree d. Then, p{x) is quasiconvex if and only 
if it can be written as 

p{x) = hi^x), (28) 
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for some nonzero ^ G M", and for some monotonic univariate polynomial h{t) of degree d. If, in 
addition, we require the nonzero component of ^ with the smallest index to he equal to unity, then 
^ and h{t) are uniquely determined by p{x). 

Proof. It is easy to see that any polynomial that can be written in the above form is quasicon- 
vex. In order to prove the converse, let us assume that p{x) is quasiconvex. By the definition of 
quasiconvexity, the closed set S{a) = {x \ p{x) < a} is convex. On the other hand, Lemma 
states that the closure of the complement of S{a) is also convex. It is not hard to verify that, as a 
consequence of these two properties, the set S{a) must be a halfspace. Thus, for any given a, the 
sublevel set S{a) can be written as {x \ ^{a^x < c{a)} for some ^(a) € M"' and c(a) G M. This of 
course implies that the level sets {x \ p{x) = a} are hyperplanes of the form {x \ £,{a)'^x = c{a)}. 

We note that the sublevel sets are necessarily nested: if a < /3, then S{a) C 5(/3). An 
elementary consequence of this property is that the hyperplanes must be collinear, i.e., that the 
vectors ^(a) must be positive multiples of each other. Thus, by suitably scaling the coefficients c(a), 
we can assume, without loss of generality, that ^(a) = ^, for some ^ G M", and for all a. We then 
have that {x \ p{x) = a} = {x \ £,^x = c{a)}. Clearly, there is a one-to-one correspondence between 
a and c(a), and therefore the value of p{x) is completely determined by S,^x. In particular, there 
exists a function h{t) such that p{x) = h{q^x). Since p{x) is a polynomial of degree d, it follows 
that h{t) is a univariate polynomial of degree d. Finally, we observe that if h{t) is not monotonic, 
then p{x) is not quasiconvex. This proves that a representation of the desired form exists. Note 
that by suitably scaling ^, we can also impose the condition that the nonzero component of ^ with 
the smallest index is equal to one. 

Suppose that now that p{x) can also be represented in the form p{x) = h{^x) for some other 
polynomial h{t) and vector ^. Then, the gradient vector of p{x) must be proportional to both ^ and 
^. The vectors ^ and ^ are therefore collinear. Once we impose the requirement that the nonzero 
component of | with the smallest index is equal to one, we obtain that ^ = ^ and, consequently, 
h = h. This establishes the claimed uniqueness of the representation. □ 

Remark. It is not hard to see that iip{x) is homogeneous and quasiconvex, then one can additionally 
conclude that h[t) can be taken to be h{t) = f^, where d is the degree of p{x). 

Theorem 4.7. For any fixed odd degree d, the quasiconvexity of polynomials of degree d can he 
checked in polynomial time. 

Proof. The algorithm consists of attempting to build a representation of p{x) of the form given in 
Proposition 14.61 The polynomial p{x) is quasiconvex if and only if the attempt is successful. 

Let us proceed under the assumption that p{x) is quasiconvex. We differentiate p{x) symboli- 
cally to obtain its gradient vector. Since a representation of the form given in Proposition 14. 61 exists . 
the gradient is of the form SIp{x) = ^h'{^x)^ where h'{t) is the derivative of h{t). In particular, the 
different components of the gradient are polynomials that are proportional to each other. (If they 
are not proportional, we conclude that p{x) is not quasiconvex, and the algorithm terminates.) By 
considering the ratios between different components, we can identify the vector ^, up to a scaling 
factor. By imposing the additional requirement that the nonzero component of ^ with the smallest 
index is equal to one, we can identify ^ uniquely. 

We now proceed to identify the polynomial h{t). For A: = 1, . . . , (i+ 1, we evaluate p{k(,), which 
must be equal to h{^'^£,k). We thus obtain the values of h{t) at d + 1 distinct points, from which 
h{t) is completely determined. We then verify that h{^'^x) is indeed equal to p{x). This is easily 
done, in polynomial time, by writing out the 0{n'^) coefficients of these two polynomials in x and 
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verifying that they are equal. (If they are not all equal, we conclude that p{x) is not quasiconvex, 
and the algorithm terminates.) 

Finally, we test whether the above constructed univariate polynomial h is monotonic, i.e., 
whether its derivative h'(t) is either nonnegative or nonpositive. This can be accomplished, e.g., 
by quantifier elimination or by other well-known algebraic techniques for counting the number and 
the multiplicity of real roots of univariate polynomials; see [5]. Note that this requires only a 
constant number of arithmetic operations since the degree d is fixed. If h fails this test, then p(x) 
is not quasiconvex. Otherwise, our attempt has been successful and we decide that p{x) is indeed 
quasiconvex. □ 

4.2.2 Pseudoconvexity of polynomials of odd degree 

In analogy to Proposition 14.61 we present next a characterization of odd degree pseudoconvex 
polynomials, which gives rise to a polynomial time algorithm for checking this property. 

Corollary 4.8. Let p{x) be a polynomial of odd degree d. Then, p{x) is pseudoconvex if and only 
if p{x) can be written in the form 



for some ^ € M" and some univariate polynomial h of degree d such that its derivative h'{t) has no 
real roots. 

Remark. Observe that polynomials h with h' having no real roots comprise a subset of the set of 
monotonic polynomials. 

Proof. Suppose that p{x) is pseudoconvex. Since a pseudoconvex polynomial is quasiconvex, it 
admits a representation h{^'^x) where h is monotonic. If h'{t) = for some t, then picking a = t ■ 
^/ll^llii '^^ have that Vp(a) = 0, so that by pseudoconvexity, p{x) is minimized at a. This, however, 
is impossible since an odd degree polynomial is never bounded below. Conversely, suppose p{x) 
can be represented as in Eq. ()29p . Fix some x, y, and define the polynomial u{t) = p[x + t{y — x)). 
Since u{t) = h{^'^x + t£/^{y — x)), we have that either (i) u{t) is constant, or (ii) u'{t) has no real 
roots. Now if Vp{x){y — x) > 0, then u'{Q) > 0. Regardless of whether (i) or (ii) holds, this implies 
that u'{t) > everywhere, so that n(l) > n(0) or p{y) > p{x). □ 

Corollary 4.9. For any fixed odd degree d, the pseudoconvexity of polynomials of degree d can be 
checked in polynomial time. 

Proof. This is a simple modification of our algorithm for testing quasiconvexity (Theorem l4.7p . The 
first step of the algorithm is in fact identical: once we impose the additional requirement that the 
nonzero component of with the smallest index should be equal to one, we can uniquely determine 
the vector ^ and the coefficients of the univariate polynomial h{t) that satisfy Eq. (j29p . (If we 
fail, p{x) is not quasiconvex and hence also not pseudoconvex.) Once we have /i(t), we can check 
whether h'{t) has no real roots e.g. by computing the signature of the Hermite form of h'{t); see [5]. 



Remark 4.1. Homogeneous polynomials of odd degree d>2> are never pseudoconvex. The reason 
is that the gradient of these polynomials vanishes at the origin, but yet the origin is not a global 
minimum since odd degree polynomials are unbounded below. 




(29) 



□ 
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4.3 Degrees that are hard 

The main result of this section is the foUowing theorem. 

Theorem 4.10. It is NP-hard to check quasiconvexity/pseudoconvexity of degree four polynomials. 
This is true even when the polynomials are restricted to be homogeneous. 

In view of Theorem 12.11 which estabhshed NP-hardness of deciding convexity of homogeneous 
quartic polynomials, Theorem 14.101 follows immediately from the following resultEl 

Theorem 4.11. For a homogeneous polynomial p{x) of even degree d, the notions of convexity, 
pseudoconvexity, and quasiconvexity are all equivalent^ 

We start the proof of this theorem by first proving an easy lemma. 

Lemma 4.12. Let p{x) be a quasiconvex homogeneous polynomial of even degree d>2. Then p{x) 
is nonnegative. 

Proof. Suppose, to derive a contradiction, that there exist some e > and x G M" such that 
p{x) = —e. Then by homogeneity of even degree we must have p{—x) = p{x) = — e. On the other 
hand, homogeneity of p implies that p(0) = 0. Since the origin is on the line between x and — x, 
this shows that the sublevel set <S(— e) is not convex, contradicting the quasiconvexity of p. □ 

Proof of Theorem We show that a quasiconvex homogeneous polynomial of even degree is 
convex. In view of implication (j27p . this proves the theorem. 

Suppose that p{x) is a quasiconvex polynomial. Define iS = {x G M" | p{x) < 1}. By 
homogeneity, for any a S M" with p{a) > 0, we have that 

e S. 



p{ay/'^ 



By quasiconvexity, this implies that for any a, b with p{a),p{b) > 0, any point on the line connecting 
a/p{ay/'^ and b/p{by/'^ is in S. In particular, consider 



a + b 



Because c can be written as 




\p{ayM + p{by/d J \p{ay/d^ 

we have that c € 5, i.e., p{c) < 1. By homogeneity, this inequality can be restated as 

p{a + b) < {p{aYl'' +p{bYl'')\ 

and therefore 

+ < ( p{aY/'^+p{bY'^ \ ^ p{a)+p{b) ^^^^ 



*A slight variant of Theorem 14 . 1 1 1 has appeared in [3]. 

^The result is more generally true for differentiable functions that are homogeneous of even degree. Also, the 
requirements of homogeneity and having an even degree both need to be present. Indeed, and x* — ix^ + 18a;^ are 
both quasiconvex but not convex, the first being homogeneous of odd degree and the second being nonhomogeneous 
of even degree. 
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where the last inequahty is due to the convexity of x"^. 

Finally, note that for any polynomial p, the set {x \ p{x) ^ 0} is dense in M" (here we again 
appeal to the fact that the only polynomial that is zero on a ball of positive radius is the zero 
polynomial); and since p is nonnegative due to Lemma 14.121 the set {x \ p{x) > 0} is dense in R". 
Using the continuity of p, it follows that Eq. ([30]) holds not only when a,b satisfy p{a),p{b) > 0, 
but for all a, b. Appealing to the continuity of p again, we see that for all a, b, p{Xa + (1 — X)b) < 
Xp{a) + (1 — X)p{b), for all A G [0, 1]. This establishes that p is convex. 

□ 

Quasiconvexity/pseudoconvexity of polynomials of even degree larger than four. 

Corollary 4.13. It is NP-hard to decide quasiconvexity of polynomials of any fixed even degree 
d > 4. 

Proof. We have already proved the result for d = 4. To establish the result for even degree 
d > 6, recall that we have established NP-hardness of deciding convexity of homogeneous quartic 
polynomials. Given such a quartic form p{x) := p{xi, . . . , Xn), consider the polynomial 

q{xi, Xn+i) = . . . , x„) + x^+i. (31) 

We claim that q is quasiconvex if and only if p is convex. Indeed, if p is convex, then obviously 
so is q, and therefore q is quasiconvex. Conversely, if p is not convex, then by Theorem 14. 11|, it is 
not quasiconvex. So, there exist points a,b,c G M^, with c on the line connecting a and b, such 
that p{a) < 1, p{b) < 1, but p{c) > 1. Considering points (a, 0), (&, 0), (c, 0), we see that q is not 
quasiconvex. It follows that it is NP-hard to decide quasiconvexity of polynomials of even degree 
four or larger. □ 

Corollary 4.14. It is NP-hard to decide pseudoconvexity of polynomials of any fixed even degree 
d > 4. 

Proof. The proof is almost identical to the proof of Corollarv 14.131 Let q be defined as in (j3ip . If 
p is convex, then q is convex and hence also pseudoconvex. If p is not convex, we showed that q is 
not quasiconvex and hence also not pseudoconvex. □ 

5 Summary and conclusions 

We studied the computational complexity of testing convexity and some of its variants, for polyno- 
mial functions. The notions that we considered and the implications among them are summarized 
below: 

strong convexity =^ strict convexity =^ convexity =^ pseudoconvexity =^ quasiconvexity. 

Our complexity results as a function of the degree of the polynomial are listed in Table 1. 

We gave polynomial time algorithms for checking pseudoconvexity and quasiconvexity of odd 
degree polynomials that can be useful in many applications. Our negative results, on the other 
hand, imply (under P^^NP) the impossibility of a polynomial time (or even pseudo-polynomial time) 
algorithm for testing any of the properties listed in Table 1 for polynomials of even degree four 
or larger. Although the implications of convexity are very significant in optimization theory, our 
results suggest that unless additional structure is present, ensuring the mere presence of convexity 
is likely an intractable task. It is therefore natural to wonder whether there are other properties of 
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property vs. degree 


1 


2 


odd > 3 


even > 4 


strong convexity 


no 


P 


no 


strongly NP-hard 


strict convexity 


no 


P 


no 


strongly NP-hard 


convexity 


yes 


P 


no 


strongly NP-hard 


pseudoconvexity 


yes 


P 


P 


strongly NP-hard 


quasiconvexity 


yes 


P 


P 


strongly NP-hard 



Table 1: Summary of our complexity results. A yes (no) entry means that the question is trivial 
for that particular entry because the answer is always yes (no) independent of the input. By P, we 
mean that the problem can be solved in polynomial time. 

optimization problems that share some of the attractive consequences of convexity, but are easier 
to check. 

Of course, NP-hardness of a problem does not stop us from studying it, but on the contrary, 
stresses the need for finding good approximation algorithms that can deal with a large number 
of instances efficiently. As an example, semidefinite programming based relaxations relying on 
algebraic concepts such as sum of squares decomposition of polynomials currently seem to be very 
promising techniques for recognizing convexity of polynomials and basic semialgebraic sets. It would 
be useful to identify special cases where these relaxations are exact or give theoretical bounds on 
their performance guarantees. 
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